Energy spectra of cosmic-ray nuclei at high energies 



H. S. Ahn 1 , P. Allison 2 , M. G. Bagliesi 3 ,L. Barbier 4 , J.J. Beatty 2 ,G. Bigongiari 3 ,T. J. 
Brandt 2 ,J.T. Childers 5 , N. B. Conklin 6 , S. Coutu 6 , M. A. DuVernois 5 , O. Ganel 1 , J. H. 
Han 1 , J. A. Jeon 7 , K. C. Kim 1 M. H. Lee 1 , P. Maestro 3 -*, A. Malinine 1 , PS. Marrocchesi 3 , 
S. Minnick 8 , S.I. Mognet 6 , S.W. Nam 7 , S. Nutter 9 , I. H. Park 7 , N. H. Park 7 , E. S. Seo 1 ' 10 , 
R. Sina 1 , P. Walpole 1 , J. Wu 1 , J. Yang 7 , Y. S. Yoon 1 - 10 , R. Zei 3 , S. Y. Zinn 1 

ABSTRACT 

We present new measurements of the energy spectra of cosmic-ray (CR) nuclei 
from the second flight of the balloon-borne experiment Cosmic Ray Energetics 
And Mass (CREAM). The instrument included different particle detectors to pro- 
vide redundant charge identification and measure the energy of CRs up to several 
hundred TeV. The measured individual energy spectra of C, 0, Ne, Mg, Si, and 
Fe are presented up to ~ 10 14 eV. The spectral shape looks nearly the same for 
these primary elements and it can be fitted to an £'- 2 - 66 ± - 04 power law in energy. 
Moreover, a new measurement of the absolute intensity of nitrogen in the 100-800 
GeV/n energy range with smaller errors than previous observations, clearly indi- 
cates a hardening of the spectrum at high energy. The relative abundance of N/O 
at the top of the atmosphere is measured to be 0.080 ± 0.025 (stat.) ±0.025 (sys.) 
at ~800 GeV/n, in good agreement with a recent result from the first CREAM 
flight. 
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Introduction 



Experimental studies of charged cosmic rays (CRs) are focused on the understanding 
of the acceleration mechanism of high-energy CRs, identification of their sources, and clar- 
ification of their interactions with the interstellar medium (ISM). The origin of CRs is still 
under debate, although direct measurements with instruments on stratospheric balloons or 
in space have provided information on their elemental composition and energy spectra, and 
indirect detection by ground-based experiments has extended the measurements up to the 
end of the observed all-particle spectrum. 

It is generally accepted that CRs are accelerated in blast waves of supernova remnants 
(SNRs), which are the only g alactic cand idates known with sufficient energy output to sus- 
tain the CR flux (jBelllll978l : [HillasJ |2005|) . Recent observation of emission of TeV gamma- 
rays from SNR RX J1713. 7-3946 by the HESS Cherenkov telescope array proved that high- 
energy charged particl es are accelerated in SNR shocks up to energies beyond 100 TeV 
( lAharonian et al.l 120041 ) . This result is consistent with expectations of the class of theoreti- 
cal models that predict the existence of a rigidity-dependent limit, above which the diffusive 
shock acceleration becomes inefficient. The maximum energy attainable by a nucleus of 
charge Z may range from ZxlO 14 eV to ~ ZxlO 17 eV depending on the model and t ypes 
of supernovae considered (ILagage fc Cesarskyill983l : lBerezhkolll996l : iParizot et al.ll2004l ). In 
this scenario, the energy spectra of elements exhibit a Z-dependent cutoff. As a consequence, 
the CR elemental composition is expected to change as a function of energy, marked by a 
depletion of low-Z nuclei in the several hundred TeV region. This mechanism has been 
proposed as a possible explanation of the steepening ("knee") in the CR energy spectrum 
(described by a power law: dN/dE oc E~~<) with a change of the spectral index from 7 ?»2.7 
to 7 ~ 3.1 observed at energies around 4 PeV. An alternative approach is adopted by models 
that relate the knee to leakage of CRs from the Galaxy. In this case, the knee is expected to 
occur at lower energies for light elements as compared to heavy nuclei, due to the rigidity- 
dependence of t he Larmor radius of CR particles propagating in the galactic magnetic field 
jHorandellbooi ). 

A detailed understanding of the propagation of CRs during their wanderings inside the 
Galaxy and their interactions with the ISM is needed to infer the injection spectra of indi- 
vidual elements at the source from the observed spectra at Earth. CR nuclei can be divided 
into primaries (i.e., H, He, C, O, Fe, etc.) that come from CR sources, and secondaries (such 
as Li, Be, B, Sc, Ti, V, etc.), which are products of the interactions of primary nuclei with 
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the ISM. The amount of material traversed by CRs between injection and observation can be 
derived from the measured ratio of secondary-to-primary nuclei, such as the boron-to-carbon 
ratio (B/C) or the ratio of sub- iron elements to iron. Likewise, their confinement time in the 
Galaxy can be determined through measurements of long-lived radioactive secondary nuclei 
( iLongairl 1 1994 ; lYanasak et al.l 1200 ll ) . Observat ions from space-bas ed experi ments like CRN 



on th e Spacelab2 mission of the Space Shuttl e (jSwordv et al.lll990f) and C2 (lEngelman et al. 
199Clh and Heavy Nuclei Experiment (HNE) fainns et alill988h onboard the HE AO 3 satel- 



lite have shown that the energy spectra of the secondary nuclei are steeper than those of the 
primaries. The escape pathlength of CRs from the Galaxy depends on the magnetic rigidity 
as R~ s , with the parameter 5 ~ 0.6. This implies that the spectra observed at Earth are 
steeper than the injection spectra, i.e., the spectral index at the source is smaller by the 
value 5. 

Those pioneering measurements were statistically limited at energies of order ~ 10 11 eV. 
Accurate direct measurements of the energy spectra of individual elements into the knee re- 
gion are needed to discriminate between different astrophysical models proposed to explain 
the acceleration and propagation mechanisms. Recently, a new generation of balloon-borne 



(Panov et al. 


2006; 


Ave et al. 


2008) 



Energetics And Mass (CREAM) experiment was designed to measure directly the elemental 
composition and the energy spectra of CRs f rom hydrogen to iron over the energy range 10 11 - 
10 15 eV in a series of flights (ISeo et al.l 12004 ). Since 2004, four instruments were successfully 



flown on long-duration balloons in Antarctica. The instrument configurations varied slightly 
in each mission due to various detector upgrades. The first flight produced important results 
extending the measurement of the relative abundances of CR secondary nuclei B and N close 
to the highest energies (~ 1.5 TeV/n) allowed by the irreducible background generated by 
the residual atmospheric overburden at balloon altitudes. The data clearly indicate that the 
escape pathlength of CRs from t he Galaxy has a power- law rigidity-dependence, with the 



parameter 5 in the range 0.5-0.6 (lAhn et al.l 120081 ) 



In this work, we present new measurements of the energy spectra of the major primary 
CR nuclei from carbon to iron, up to a few TeV/n made with the second CREAM flight 
(CREAM-II). A new measurement of the nitrogen intensity with unprecedented statistics is 
also presented up to 800 GeV/n. The measured N/O ratio confirms the results of CREAM-I. 
The procedure used to analyze the data and reconstruct the CR energy spectra is described 
in detail together with an assessment of the systematic uncertainties. 
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2. The CREAM-II instrument 

The instrument for the second flight included: a redundant system for particle identi- 
fication, consisting (from top to bottom) of a timing-charge detector (TCD), a Cherenkov 
detector (CD), a pixelated silicon charge detector (SCD), and a sampling imaging calorime- 
ter (CAL) designed to provide a measurement of the energy of CR nuclei in the multi-TeV 
region. Figure [1] depicts the detector arrangement. 

The TCD is made of two orthogonal layers of four 5 mm-thick plastic scintillator paddles 
each, covering an area of 120x120 cm 2 . The paddles are read out by fast photomultiplier 
tubes (PMTs) via twisted-strip adiabatic light guides. A custom design of the electronics 
combines fast peak detection with a threshold timing measurement of the PMTs signals in 
order to determine each element's charge with a resolution < 0.35 e (in units of the electron 
charge) and to discriminate against albedo particles. A hodoscope (S3) of 2x2 mm 2 square 
scintillating fibers read out by two PMTs is located just above the calorimeter, and provides 
a reference time for TCD timing rea d out. Further details about the TCD design and per- 



formance can be found in lAhn et al.l (120091 ) . The CD is a 1 cm-thick plastic radiator with 
1 m 2 surface area, instrumented with eight PMTs viewing wavelength shifting bars placed 
along the radiator edges. It is used to flag relativistic particles and to provide a charge 
determination complementary to that of the TCD. 

The SCD is comprised of two layers of 156 silicon sensors each. Each 380 /zm-thick sen- 
sor is segmented into an array of 4x4 pixels, each of which has an active area of 2.12 cm 2 



( IPark et al.l 120071 ). The sensors are slightly tilted and overlap each other in both lateral 
directions, providing a full coverage in a single layer, with a 77.9 x 79.5 cm 2 area. The 
dual layers, about 4 cm apart, cover an effective inner area of 0.52 m 2 with no insensitive 
regions. The overall SCD vertical dimension, including the lad ders, the mechan ical support 



structure and the electromagnetic shielding cover, is 97.5 mm (INam et al.l 120071 ) . The read- 
out electronics of the 4992 pixels is based on a 16-channel ASIC chip followed by 16-bit 
analog-to-digital converters (ADCs), providing a fine charge resolution over a wide dynamic 
range from hydrogen to nickel. 

The CAL is a stack of 20 tungsten plates, each 1 radiation length (X ) thick and with surface 
area 50x50 cm 2 , interleaved with active layers of 1 cm- wide and 50 cm-long scintillating-fiber 
ribbons. Each ribbon is built by gluing together 19 scintillating fibers of 0.5 mm diameter 
each. The light signal from each ribbon is collected by means of an acrylic light mixer coupled 
to a bundle of 48 clear fibers. This is split into three sub-bundles (with 42, 5 and 1 fibers, 
respectively), each feeding a pixel of a hybrid photodiode (HPD). This optical division is 
used to match the wide dynamic range of the calorimeter to that of the front-end electronics, 
providing three readout scales (low, medium, and high) with different sensitivity. A total of 
2560 channels are readout from 40 HPDs powered in groups of 5 units. 
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The longitudinal and lateral segmentations of the CAL correspond to 1 Xq and 1 Moliere 
radius, respectively. The ribbon planes are alternately oriented along orthogonal directions 
to image in three dimensions, the development of showers produced by the incoming nu- 
clei which interact inela stically in the ~ 0.46 Aj„ r thick carbon target preceding the CAL 
( jMarrocchesi et al.ll2005l ). The finely grained sampling calorimeter can track, by reconstruct- 
ing the shower axis, the incident particle trajectory with enough accuracy to determine which 
segment of each charge detector it traversed, and thereby allowing discrimination against seg- 
ments hit by the backscattered particles produced in the shower. This feature is essential 
for a reliable and unambiguous charge identification. 



3. Balloon flight and instrument performance 

The second CREAM instrument was launched from McMurdo station in Antarctica on 
2005 December 16th. The balloon circled the South Pole twice at a float altitude between 
35 and 40 km. The flight was terminated on 2006 January 13th, 28 days after the launch. 
During the whole flight, the instrument housekeeping system monitored online temperatures, 
currents, voltages of the front-end electronics boards, and the high-voltage supply values of 
the HPDs and PMTs. The performance of all sub-detectors met the design technical speci- 
fications. The CAL and TCD high-voltage systems worked successfully at low atmospheric 
pressure. The thermal behavior of the various detectors was found to be in good agreement 
with expectations from thermal models. The temperature of the instrument crates stayed 
within the required operational range with daily variation of a few °C, depending on the 
inclination of the Sun during each 24 hr cycle. Front-end electronics channels showed stable 
gains over the whole flight. Pedestal values of the 2560 CAL and 4992 SCD channels were 
collected automatically every 5 minutes in order to monitor their drift as a function of the 
detector temperature and to perform accurate pedestal subtraction. 1.7% of the SCD and 
0.7% of the CAL channels were noisy or inefficient and were masked off. 
During the flight, the data acquisition was enabled whenever a shower developing in at least 
six planes was detected in the CAL (CAL trigger) and/or a relativistic CR with Z > 2 was 
identified by the TCD and CD (TCD trigger). The high-energy data, i.e., the CAL triggered 
events, were transmitted via Tracking and Data Relay Satellite System (TDRSS), while the 
low-energy data, i.e., events triggered only by the TCD and not by the CAL, were recorded 
on an onboard disk. A total of 57 GB of data were collected. In order not to saturate the 
disk space too early in the flight, the highly abundant low-energy particles triggered only by 
the TCD were prescaled, so that only a fraction of these events were recorded. This fraction 
of recorded TCD triggers (prescale factor) could be set by the data acquisition program and 
changed during the flight depending on the trigger rate. On average, only one every six TCD 
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triggered events was retained. A detailed description of the instrument performance during 



the second flight may be found in iMarrocchesi et al.l (120081 ). 



4. Data analysis 

The present analysis uses a subset of data collected from December 19th to January 
12th, under stable instrument conditions, representing a total of 24 days' worth of data 
taking (e.g., removing periods early in the flight and at the very end, when the instrument 
was either being adjusted and optimized, or prepared for flight termination). The first step 
in the analysis is to correct the raw data for the gain variations among the channels and 
the drift of the pedestals with temperature. Then the trajectory of each CR through the 
instrument is reconstructed and each event assigned a charge and an energy. At this point, 
the events of each charge species are sorted into energy intervals and the number of counts 
in each bin is corrected for several effects in order to compute the differential intensity. The 
various steps of the analysis procedure are described in detail in the following subsections. 



4.1. Trajectory reconstruction 

To determine the arrival direction of a CR, the axis of the shower imaged by the calorime- 
ter is measured with the following method. In each CAL plane, an offline software algorithm 
searches for clusters of adjacent cells with pulse heights larger than a threshold value, cor- 
responding to an energy deposit of about 10 MeV. The cluster formed by the cell with the 
maximum pulse height and its two neighbors defines a candidate track point. Its coordinates 
are computed as the center of gravity of the cluster. A track is formed by matching at least 
three candidate track points sampled in each CAL view. The shower axis parameters are 
calculated by a linear x 2 fit of the track. Track quality is assured by requiring a value of 
X 2 < 10. The reconstructed shower axis is projected back to the top of the instrument in or- 
der to determine which SCD pixels and TCD paddles were hit by the incoming CR. Particles 
showering in the CAL, but crossing neither the TCD scintillator paddles nor the SCD planes 
are rejected. Otherwise, if the particle crosses the SCD, a circle of confusion with a 3 cm 
radius centered on the impact point of the track is traced in each SCD layer. The value of 
the radius is equal to about three times the estimated uncertainty on the impact point. The 
hit pixels within each circle of confusion are scanned and the pixel with the maximum pulse 
height is selected for the charge assignment procedure. In order to improve the trajectory 
determination and the accuracy of the pathlength correction, the two selected pixels, one 
per SCD layer, are added to the candidate points and the track is fitted again. In this way, 
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the intersection of the particle trajectory with the top SCD layer is measured with a spatial 
resolution of 7 mm rms. The incidence direction is determined with an accuracy of better 
than 1° as estimated with Monte Carlo (MC) simulations. 

This procedure reconstructs well the trajectories of over 95% of the CR particles with energy 
> 3 TeV triggered by the CAL. It can also identify the TCD triggered events which pro- 
duce a shower in the CAL, but with an energy release insufficient to satisfy the CAL trigger 
condition. 



4.2. Charge assignment 

The identification of the charge Z of a CR relies on two independent samples of its 
specific ionization dE/dx provided by the SCD layers. Since dE/dx depends on Z 2 and for 
solid absorbers it is nearly constant in the relativistic regime, the particle charge Z can be 
determined by measuring the amount of ionization charge a CR produces when traversing 
the silicon sensors. In order to avoid charge misidentification of the incident nucleus, caused 
by the backscattered shower particles reaching the SCD, and to obtain high-purity samples 
of CR elements, the signals, S top and Sbottom, of the two pixels selected by the tracking 
algorithm are compared and required to be consistent within 20%. If this coherence cut 
is satisfied, they are corrected for the pathlength (estimated from the track parameters) 
traversed by the particle in the silicon sensors, which is proportional to 1/ cos (9), 9 being 
the angle between the particle trajectory and the axis of the instrument. A reconstructed 
charge Z rec is assigned to the particle by combining the ionization signals matched with the 
track on both layers, according to the formula: 



7 i Stop + Sbottom COS (9) , . 

rec ~\i 2 X — ^ — ^ ' 

where sq is a calibration constant, inferred from the flight data, to convert the pulse heights 
into an absolute charge scale. The reconstructed charge distribution is shown in Figure EJ 
The elemental range from boron to silicon is fitted with a multi-Gaussian function (Figure 



2(a)). From the fit, a charge resolution a is estimated as: 0.2 e for C, N, O and ~ 0.23 e 
for Ne, Mg, Si. Though a multi-Gaussian fit in the elemental region from sulfur to nickel 
cannot be done due the limited statistics, the charge resolution for iron is estimated to be 
~ 0.5 e from the width of the Fe peak. A 2er cut around the mean charge value is applied 
to select samples of C, O, Ne, Mg and Si. A cut of Z rec ± la is imposed for Fe due to the 
corresponding lower SCD resolution, while for N a cut of Z rec ±1.5er and 10% coherence level 
between the selected pixels signals, are used to avoid contamination from more abundant 
adjacent charges. 
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4.3. Energy measurement 

Due to the limitations on mass and size imposed to the payload, a balloon-borne ex- 
periment cannot employ a total containment hadronic calorimeter to measure the energy of 
CR nuclei. An alternative and workable technique is to use a thin ionization calorimeter 
to sample the electromagnetic core of the hadronic cascade initiated by a CR interacting 
in a target preceding the calorimeter. In fact, though a significant part of the hadronic 
cascade energy leaks out of the calorimeter, the energy deposited in the calorimeter by the 
shower core still scales linearly with the incident particle energy, albeit with large event-to- 
event fluctuations. As a result, the energy resolution is poor by the standards of hadron 
calorimetry in experiments at accelerators. Nevertheless, it is sufficient to reconstruct the 
steep energy spectra of CR nuclei with a nearly energy independent resolution. 
The CREAM calorimeter was designed to measure CRs over a wide energy range, from tens 
of GeV up to about 1 PeV. It is characterized by a very small sampling fraction; only about 
0.1% of the CR energy is converted into visible signals in the CAL active planes and collected 
by the photosensors. Three different gain scales (low, medium and high) were implemented 
for each ribbon in order to prevent saturation of the readout electronics and thereby to en- 
sure a linear detector response up to very high energies. 

The total energy deposited in the calorimeter by an interacting nucleus, is measured by sum- 
ming up the corrected pulse heights of all the cells. The cells are equalized for non-uniformity 
in light output and gain differences among the photodetectors using a set of calibration con- 
stants from accelerator beam tests (described below), scaled according to the high- voltage 
settings during the flight. Low and medium gain scales are inter-calibrated using the flight 
data; the medium-range signal of each ribbon is used whenever the corresponding low-range 
signal deviates from linearity due to saturation. 

The CREAM calorimeter was tested, equalized, and calibrated pre-flight at the CERN 



SPS accelerator with beams of electrons and protons up to a few hundred GeV (lAhn et al. 



20041 ) . Its response to relativistic nuclei was also studied by exposing the calorimeter to 
nuclear fragments from a 158 GeV/n primary indium beam. The detector was found to have 
a linear response up to the maximum beam energy, equivalent to about 8.2 TeV particle 
energy, and a nearly flat reso lution (around 30%), at energies above 1 TeV for all heavy nuclei 



with Z>5 (lAhn et al.l 120061 ) . Above the maximum beam energy, the calorimeter response 
was extrapolated using MC simulations, which predict a linear behavior up to hundreds of 
TeV. 
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4.4. Energy deconvolution 

Once each CR is assigned a charge and energy, the reconstructed particles of each 
nuclear species are sorted into energy intervals commensurate with the rms resolution of the 
calorimeter. Due to the finite energy resolution of the detector, the measured number of 
events in each energy bin must be corrected for overlap with the neighboring bins. This 
unfolding procedure requires solving a set of linear equations 

n 

M % = Y^^jNj i = l,...,m (2) 

3=1 

relating the "true" counts Nj in one of n incident energy bins to the measured counts Mj 
in one of m deposited energy bins. A generic element of the mixing matrix Ay represents 
the probability that a CR particle, carrying an energy corresponding to a given energy bin 
j, produces an energy deposit in the calorimeter falling into bin i inste ad. A detailed MC 



simulation of the instrument, based on the FLUKA 2006.3b package (IFasso et al.ll2005l ). 
was developed to estimate the unfolding matrix. Sets of nuclei, generated isotropically 
and with energies chosen according to a power-law spectrum, are analyzed with the same 
procedure used for the flight data. Each matrix element Ajj is calculated by correlating 
the generated s pectrum with the distribution of the deposited energy in the calorimeter 



(IZei et al.l 120081 ) . In order to get a reliable set of values of the unfolding matrix for each 
nucleus, the MC simulation is finely tuned to reproduce both flight data and the calibration 
data collected with accelerated particle beams. The agreement of the MC description with 
the real instrument behavior was carefully checked. As an example, in Figure [3] the response 
of the calorimeter to carbon nuclei from the flight data is compared with an equivalent set 
of simulated events. 

The elements of the unfolding matrix depend on the spectral index assumed for the generatio n 



of MC events. For each nucleus we used the spectral index value reported by lHorandell (120031 ). 
obtained by combining direct and indirect CR observations. In order to check the stability 
of our final results, the unfolding procedure was repeated by scanning several values of the 
spectral index in the range ±0.2 around the reference value. No significant bias of the final 
results was observed. 



4.5. Corrections for interactions in the atmosphere and the detector 

CRs may undergo spallation reactions when traveling through the atmosphere above 
the balloon or crossing the instrument material of CREAM above the SCD. This results in 
the generation of secondary nuclei with lower charge than the parent's, and of about the 
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same energy per nucleoli. Then, for each nuclear species, inelastic interactions lead to a 
loss of primary particles but also a gain of secondaries produced by heavier nuclei. For a 
correct determination of the CR fluxes, the data must be corrected for this effect. A two-step 
procedure is used which takes into account separately the CR pathlength in air (correction 
to the top of the atmosphere (TOA)) and the amount of instrument material encountered 
by the particle before reaching the SCD. The latter leads to a correction to the top of the 
instrument (TOI). For each step the corrections are calculated by solving the matrix equation 
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which describes the transport process of the nuclei from an incident level X (e.g., the top of 
the instrument) to a final level T (e.g., the SCD). The vector Nj represents the number of 
nuclei of each species counted in a given energy interval, while the vector Njr is the number 
of counts corrected for the charge- changing interactions in the materials between the levels 
X and T. The diagonal elements of the transport matrix represent the survival probabilities 
of each species Z in the passage from X to J 7 , while a generic off-diagonal element Fz 2 ^z 1 
gives the fraction of nuclei Z 2 which spallate to produce nuclei of charge Z\ (< Z 2 ). 
The TOI transport matrix is obtained from the MC simulation of the instrument. The mass 
thickness of the CREAM-II materials above the SCD is estimated to be ~ 4.8 g cm~ 2 . The 
fraction of surviving nuclei, i.e., the diagonal matrix elements, ranges from 81.3% for C to 
61.9% for Fe. These MC values are compared with the ones calculated using the sem i empir - 
ical total cross-section formulae for nucleus-nucleus reactions given in ISivher et al.l (119931 ). 
Good agreement, within 1%, is found as shown in Figure HI 

The TOA matrix elements are calculated by simulating with FLUKA the atmospheric over- 
burden during the flight (3.9 g cm -2 on average). Since the probability of interaction depends 
on the nucleus pathlength in air, the simulated events are generated according to the mea- 
sured zenith angle distribution of the reconstructed nuclei, which is peaked around 20°. 
Survival probabilities ranging from 84.2% for C to 71.6% for Fe are found (Figure ED- They 
turn out to be in good agreement w i th the values ca lculated from two different empirical 
parameterizations (ISivher et al.l 1 19931 ; iKox et al.lll987l ) of the total fragmentation cross sec- 
tions of heavy relativistic nuclei in targets of different mass. 

The off-diagonal elements Fzi-+Z\ of both the TOI and TOA transport matrices are of order 
l%-2%. The comparison between values obtained from different models is discussed in Sec- 
tion O 
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4.6. Absolute differential intensity calculation 



The absolute differential intensity dN/ dE at an energy E is calculated for each elemental 
species according to the formula 



dN 
dE 



Ni 



1 



AEa ex SttxT 



(4) 



where iVj is the number of events obtained from the unfolding algorithm and corrected to the 
TOA, AEi is the energy bin width, T is the exposure time, Sfl the geometric factor of the 
instrument, and e the efficiency of the analysis selection. Since a significant fraction of the 
reconstructed events at energies <3 TeV did not satisfy the CAL trigger condition and were 
triggered only by the TCD, the number of these events in the first energy intervals must be 
further divided by the prescale factor (1/6 on average during the flight) of the TCD-based 
trigger system. 



Median energy. According to the work of iLafferty fc Wyattl (119951 ). each data point, 
i.e., each differential intensity value, is centered at a median energy E, defined as the energy 
at which the measured spectrum is equal to the expectation value of the "true" spectrum. 
For a power-law spectrum (E~^) with spectral index 7, the median energy is calculated as 



E 



TTll-7 

^2 



E[ 



1-7 



(E 2 -E 1 ) (I-7) 



(5) 



where Ei and E 2 are the lower and upper limits of a given energy bin. If the highest energy 

1 

interval is not limited, i.e., E 2 — > 00, the median energy is E = 2~<- 1 E\. In this case, the 
energy interval is an integral bin and its width (to be used in Equation H|) is calculated as 



AE 



E x 



2— ( 7 _ 1) 



(6) 



Since 7 is not known a priori, an iterative procedure is performed to c ompute the med ian 
energy of each bin. For each element, the spectral index reported by iHorandell (120031 ) is 
used as the initial value of 7. We found that E depends weakly on 7 for a variation of ± 0.1 
around the initial value. 

Geometric factor. The geometric factor SQ is estimated from MC simulations by count- 
ing the fraction of generated particles entering the trigger-sensitive part of instrument. Two 
groups of events are distinguished as shown in Figure [TJ CRs crossing both the TCD and 
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the SCD before impinging onto the CAL (named "golden" events), and particles entering 
the instrument through the side (i.e., crossing neither the TCD nor the CD) and traversing 
only the SCD and the CAL. The estimated acceptance is (0.19 ± 0.01) m 2 sr in the first case, 
while including both event types yields a value of (0.46 ± 0.01) m 2 sr. Both values are charge 
and energy independent. For the second group of events, the corrections to the TOI are not 
necessary. 

Live Time. During the flight, the total time and the live time T (i.e., the time during 
which the data acquisition was available for triggers) were measured by means of a pair of 
48-bit counters incorporated in the housekeeping system onboard, and providing better than 
4 microseconds resolution. The selected data set amounts to a live time of 1,454,802 seconds, 
close to 75% of the total time of data taking. 

Reconstruction efficiency. The overall efficiency e, including the trajectory reconstruc- 
tion and the charge selection efficiencies, is estimated from MC simulations as a function of 
the particle energy. This efficiency curve is characterized by a steep rise below 2 TeV and 
reaches a constant value at energies >3 TeV for all nuclei. The rise is caused by the high 
sparsification threshold (~10 MeV) of the CAL cells, which limits the CAL tracking accu- 
racy at low energy and consequently the charge identification capability of the instrument. 
For C, O, Ne, Mg, and Si, the plateau value in the efficiency curve is about 80% when con- 
sidering only the "golden" events, and decreases to 65% when including also events at large 
angle within the SCD-CAL acceptance. Instead, the more conservative charge cuts applied 
to select N and Fe samples (see Section 14.21) yield a reduced plateau efficiency of about 65% 
and 35% (for "golden" events), respectively. 



5. Results 

The differential intensities at the TOA as measured by CREAM-II for the elements C, 
N, O, Ne, Mg, Si, and Fe are given as a function of the kinetic energy per nucleon in Table 
[[J together with the corresponding number of raw observed particles. The quoted error bars 
are due to the counting statistics only: the upper and l ower limi t s are computed as 84% 



confidence limits for Poisson distributions as described by iGehreld ( 119861 ). 
The energy spectra of the major primary CR nuclei from carbon to iron are plotted as a 
function of the kinetic energy per particle in Figure EJ To emphasize the spectral differences, 
the differential intensities are multiplied by E 2 ' 5 and plotted as a function of the kinetic 
energy per nucleon in Figure [7J The error bars shown in the figures represent the sum in 
quadrature of the statistical and systematic uncertainties. The assessment of systematics is 
discussed in detail in the following section. The particle energy range covered by CREAM-II 
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extends from around 800 GeV up to 100 TeV. The absolute intensities are presented without 
any arbitrary normalization to previous data and cover a range of six decades. 
The energy spectrum of nitrogen is shown in Figure [SJ The statistical and systematic errors 
in the differential intensity are shown separately to point out that the systematic uncer- 
tainties from the corrections for secondary particle production in the atmosphere and the 
instrument, are particularly relevant for this measurement. In practice, they impose a limit 
on the maximum energy at which measurements of nitrogen in CRs might be pursued on a 
balloon flight. 

The elemental abundances with respect to oxygen are calculated from the differential in- 
tensities at the TOA and are shown as a function of energy in Figure [9j Energy bins at 
high energy are merged together in order to reduce the statistical error in the ratio to an 
acceptable level. 



5.1. Estimate of the systematic errors 



The main systematic uncertainties in the differential intensity stem from the recon- 
struction algorithm as well as from the TOI and TOA corrections. To estimate the first, 
we scan a range of thresholds around the reference value of each selection cut and derive 
the corresponding fluxes. In particular, the analysis procedure has been repeated by varying 
the fiducial area of the SCD and CAL, the coherence level between the signals of the SCD 
pixels crossed by the CR particle, and the limits of charge selection cut for each element. 
Comparing the results with the reference, a systematic fractional error in the differential 
intensity, arising from the reconstruction procedure, is estimated to be of order 10% below 
a particle energy of 3 TeV and 5% above. 

The second source of systematic errors comes from the uncertainties in the nucleus-nucleus 
charge- changing cross sections used to calculate the instrument and atmospheric corrections. 
While there is a general agreement between different parameterizations of the total cross sec- 
tions, leading to very similar values of the survival fractions of CRs in the atmosphere and 
in the instrument (Figures H] and Ej), the partial charge-changing cross sections reported in 
the literature are measured with quoted errors of order 10%-15%, depending on the specific 
nucleus-nucleus interaction, and can differ by up to ~30%. This may result in significantly 
different values of the off-diagonal elements Fz % -*z x in the transport matrix. For instance, the 
values we calculated for the fraction of oxygen nuclei which spallate in the atmosphere pro- 
d ucing nitrogen nuc lei (F q ^n), are respectively: 1.7% from the partial cross-section model 
of iTsao et al.l ( 19931); 1-9% acco rding to the parameterization of iNilsen et al.l (119951 ); 1.5% 
from FLUKA (IFasso et al.ll2005l ). In order to estimate how these uncertainties in the correc- 
tions for secondary nuclei production in the atmosphere and the instrument affect our final 
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results, we developed the following method. Using the same sample of analyzed data for each 
element, we varied the values of the charge- changing probabilities in the range allowed by the 
different models and derived, for each set of values, the corresponding energy spectrum. By 
comparing the results with the reference spectrum, the systematic error for the uncertainty 
in correcting the differential intensity to the top of instrument is estimated to be 2% for the 
primary elements. For nitrogen, a 15% error is assigned because of the large contamination 
of O nuclei which spallate into N. Similar values are estimated for the systematic error from 
the uncertainties in the atmospheric secondary corrections. 

The systematic error in the energy scale, as measured by the CAL, derives from the accuracy 
with which the simulated CAL response corresponds to the real behavior of the calorimeter 
at energies not covered by beam test measurements. We estimate this accuracy at a 5% level, 
which corresponds, given the CAL linearity in the energy range considered in this work, to 
an energy-scale uncertainty of 5%. 



Discussion 



The CREAM-II results for primary nuclei ar e found to be in good a greeme nt with previ- 



ous m easurements of space-based ( HEAP 3 C2 ( Engelman et al.lll990T). CRN (jMiiller et al. 



19911 1) and balloon-borne (ATIC-2 flPanov et al. 



2006h. TRACER jAve et aHl2008h ) experi- 



ments (see Figures El and [7]) . All the elements appear to have the same spectral shape, which 
can be described by a single power law in energy _E~ 7 . The spectral indices fitted to our 
data (Table [2]) are very similar, indicating that the intensities of the more abundant (evenly) 
charged heavy elements have nearly the same energy dependence. Our observations, based 
on a calorimetric measurement of the CR energy, confirm the results recently reported by 
the TRACER collaboration (Figure [TUT) , using completely different techniques (Cherenkov, 
specific ionization in gases and transition radiation) to determine the particle energy. The 
weighted average of our fitted spectral indices is 7 = 2.66 ± 0.04, consistent, within error, 
with the value of 2.65 ± 0.05 obtained from a fit to the combined CRN and TRACER data 



(jAve et al.ll2008l ). The great similarity of the spectral indices suggests that the same mech- 
anism is responsible for the source acceleration of these primary heavy nuclei. 
Although all the elemental spectra can be fitted to a single power law in energy, hardening 
(flattening) of the spectra above ~ 200 GeV/n is apparent from the overall trend of the data 
(Figure UJ). A detailed analysis to investigate possible features in the spectra is in progress 
to evaluate the significance of a change in the spectral index. 

CREAM-II data extend the energy range spanned by the previous elemental abundance 
measurements up to ~ 1 TeV/n (Figure [9]). The C/O ratio m easured by CREAM-II is con- 
sistent with the one reported by CREAM-I (lAhn et al.l 120081 ) within the overlapping energy 
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region covered by the two nights and slightly higher than the CRN result. The Ne/O and 
Mg/O ratios confirm, up to 800 GeV/n, the nearly flat trend of HE AO 3 C2 low-energy 
data, while the fractions of Si and Fe with respect to O seem to increase with energy This 
could represent the clue of a possible change in composition of CRs at high energies with an 
enhancement of heavier elements, as expected from SNR shock diffusive acceleration theo- 
ries. Nevertheless, the statistics of the data are too limited for a definitive conclusion at this 
time. 

Unlike the primary heavy nuclei, nitrogen is mostly produced by spallation in the inter- 
stellar medium , but it also has a primary contribution of order 10%, as recently measured 



by CREAM-I (lAhn et al.ll2008l ). The measurement of nitrogen intensity at high energy is 
challenging because it requires an excellent charge separation between N nuclei and the much 
more abundant C and O neighbors, as well as a complete understanding of the corrections for 
secondary N particle production in the atmosphere and the instrument. The nitrogen data 
collected by CREAM-II are more statistically significant at high energy than any previous 
observation. We notice that, above 100 GeV/n, the N spectrum flattens out from the steep 
decline which characterizes the energy range 10-100 GeV/n (Figure [H]). The combined CRN 
and CREAM-II data at energies higher than 200 GeV/n, are well fitted to a power-law in 
energy with spectral index 2.61 ±0.10, consistent within the error with the average value 7 
of the primary nuclei. This supports the hypothesis of the presence of two components in 
the cosmic nitrogen flux. In fact, since the escape pathlength decreases rapidly with energy, 
the secondary nitrogen component is expected to become negligible at high energy, where 
only the primary component should survive. This might result in a change of the spectral 
slope, as observed. 

As a cross-check to the measurement reported here, we calculate the relative abundance ratio 
of N/O and compare it with the CREAM-I observation (Figure [9]), in which the energy of 
CRs was measured with a transition radiation detector. A N/O ratio (at the TOA) equal 
to 0.080 ± 0.025 (stat.) ±0.025 (sys.) is measu red at an energy of ~800 GeV/n, in good 



agreement with recent results from CREAM-I (lAhn et al.l 120081 ) 



7. Conclusions 

CREAM-II carried out measurements of high-Z CR nuclei with an excellent charge 
resolution and a reliable energy determination. The use of a dual layer of silicon sensors 
allowed us to unambiguously identify each individual element up to iron and to reduce the 
background of misidentified nuclei. We demonstrated the feasibility of measuring the particle 
energy up to hundreds of TeV by using a thin ionization sampling calorimeter, with sufficient 
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resolution to reconstruct accurately the CR spectra. 

The energy spectra of the major primary heavy nuclei from C to Fe were measured up to 
~ 10 14 eV and found to agree well with earlier direct measurements. All the spectra follow a 
power law in energy with a remarkably similar spectral index 7 = — 2.66±0.04, as is plausible 
if they have the same origin and share the same acceleration and propagation processes. 
A new measurement of the nitrogen intensity in an energy region so far experimentally 
unexplored indicates a harder spectrum than at lower energies, supporting the idea that 
nitrogen has secondary as well as primary contributions. These results provide new clues 
for understanding the CR acceleration and propagation mechanisms, but further accurate 
observations with high statistics beyond 10 14 eV are needed to finally unravel the mystery 
of CR origin. 
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Table 1. Differential Intensities Measured with CREAM-II 



Element Energy Range Kinetic Energy No. of Events Intensity 

(GeV/ra) E (GeV/n) (m 2 sr s GeV/n)- 1 



c 


66 - 


117 


86.2 


171 


(1.67 ± 0.13)xl0- 4 


Z=6 


117 - 


175 


141.6 


139 


(3.37 ± 0.29) xlO- 5 




175 - 


304 


228.2 


215 


(9.42 ± 0.64) xl0~ 6 




304 - 


529 


397.2 


134 


(2.28 ± 0.20) xl0~ 6 




529 - 


919 


691.0 


66 


(5.17 ± 0.63)xl0- 7 




919 - 


1597 


1201.5 


29 


(1.31 ± 0.24)xl0~ 7 




1597 - 


2776 


2088.7 


13 


' (3.6l^)xl0- 8 




2776 - 


4824 


3630.5 


7 


(1.04 l^)xl0- 8 




4284 - 




7415.1 


5 


(1.49 ^;l)xl0- 9 


N 


56 - 


175 


95.7 


86 


(7.28 ± 0.79) xl0~ 6 


Z=7 


175 - 


529 


295.3 


63 


(5.39 ± 0.68) xl0~ 7 




529 - 




826.3 


19 


(3.8 ±S;5)xlO" 8 





49 - 


87 


64.0 


224 


(3.48 ± 0.23) xlO" 4 


Z=8 


87- 


175 


120.8 


371 


(4.82 ± 0.25) xl0~ 5 




175 - 


304 


228.0 


241 


(7.33 ± 0.47) xl0~ 6 




304 - 


529 


396.9 


127 


(1.90 ± 0.17)xl0- 6 




529 - 


919 


690.5 


57 


(4.51 ± 0.60) xl0~ 7 




919 - 


1597 


1200.5 


25 


(1.25 ± 0.25)xl0~ 7 




1597 - 


2776 


2087.0 


10 


^ (2.9+^)xl0- 8 




2776 - 


4824 


3627.4 


6 


(8.3 ^|)xl0- 9 




4284 - 




7287.1 


3 


(8.4 ^; 7 )xl0- 10 


Ne 


40 - 


58 


47.0 


37 


(1.20 ± 0.20)xl0" 4 


Z=10 


58 - 


101 


74.8 


65 


(2.59 ± 0.32) xl0~ 5 




101 - 


175 


130.7 


67 


(5.53 ± 0.68) xl0~ 6 




175 - 


304 


227.9 


36 


(1.06 ± 0.18)xl0~ 6 




304 - 


529 


396.7 


19 


(2.86 ^)xlO- 7 




529 - 


919 


690.0 


8 


(6.5 + ^)xl0- 8 




919 - 


1597 


1199.8 


4 


(1.8 ^)xlO- 8 




1597 - 


2776 


2085.6 


2 


(3.7^)xl0- 9 




2776 - 




4150.9 


3 


(2.2 ± 2 i)xl0- 9 
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Table 1 — Continued 



Element 


Energy ] 


Range 


Kinetic Energy 


No. of Events 


Intensity 




(GeV 


In) 


E (GeV/n) 




(m 2 sr s GeV/n)- 1 




Mg 


20 - 


40 


27.0 


48 


(4.68 ± 0.67) xlO^ 4 


Z=12 


40 - 


58 


47.0 


45 


(1.07 ± 0.16)xl0- 4 




58 - 


101 


74.9 


103 


(2.91 ± 0.29) xl0~ 5 




101 - 


175 


130.8 


94 


(5.41 ± 0.56) xl0~ 6 




175 - 


304 


228.0 


51 


(1.42 ± 0.20) xl0~ 6 




304 - 


529 


397.0 


26 


(3.72 ± 0.73) xl0~ 7 




529 - 


919 


690.6 


12 


(8.5 ^)xl0~ 8 




919 - 


1597 


1200.8 


4 


(1.9 +^)xl0~ 8 




1597 - 


2776 


2087.5 


3 


(6.7^)xl0- 9 




2776 - 




4215.9 


1 


(6.0 ^ 6 8 °)xl0- 10 


Si 


20 - 


40 


27.0 


66 


(4.54 ± 0.56) xlO" 4 


z=u 


40 - 


58 


47.0 


64 


(1.35 ± 0.17)xl0" 4 




58 - 


101 


74.9 


116 


(3.03 ± 0.28) xlO" 5 




101 - 


175 


130.8 


89 


(5.04 ± 0.53) xl0~ 6 




175 - 


304 


228.0 


54 


(1.50 ± 0.20)xl0~ 6 




304 - 


529 


397.0 


24 


(3.88 ± 0.79) xlO" 7 




529 - 


919 


690.6 


10 


(1.15 1°; 48 ) xlO- 7 




919 - 


1597 


1200.7 


3 


(1.9 ±J 9 7 )xlO- 8 




1597 - 




2418.3 


2 


(3.2 ;?j)xio- 9 


Fe 


18 - 


40 


25.4 


108 


(5.19 ± 0.56) xlO" 4 


Z=2Q 


40 - 


58 


47.0 


88 


(8.38 ± 0.89) xlO" 5 




58 - 


101 


74.9 


70 


(2.06 ± 0.25) xlO" 5 




101 - 


175 


131.0 


27 


(4.46 ± 0.85) xlO" 6 




175 - 


304 


228.2 


15 


(1.36 l^)xl0- 6 




304 - 


529 


397.2 


9 


(4.0 + 1|) xlO- 7 




529 - 


919 


690.9 


4 


(1.09 l°;™)xio- 7 




919 - 




1406.9 


4 


(1.6 ^)xlO- 8 
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Table 2: Best fit values of the spectral index 7 of each element, as measured by CREAM-II, 
and x 2 /ndf of the fit 



Element 


7 


X 2 /ndf 


C 


2.61 ±0.07 


3.3/7 





2.67 ±0.07 


4.8/7 


Ne 


2.72 ±0.10 


4.2/7 


Mg 


2.66 ±0.08 


1.9/8 


Si 


2.67 ±0.08 


4.9/7 


Fe 


2.63 ±0.11 


4.2/6 
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Fig. 1. — Schematic view of the CREAM-II instrument. 
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Fig. 2. — Charge histogram obtained by the SCD in the elemental range (a) from boron to 
silicon and (b) from sulfur to nickel. The distributions are only indicative of the SCD charge 
resolution and the relative elemental abundances are not meaningful. The peaks from B to 
Si are fitted to a multi-Gaussian function. 
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Fig. 3. — Energy deposited in the calorimeter by a selected sample of carbon nuclei. Simu- 
lated (histogram) and real (dots) events are shown. 
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Fig. 4. — Survival probabilities of CR nuclei traversing the materials between the TOI and 
the upper SCD layer. For each elemental speci es, the correction is estimated in two different 
ways: (filled circles) from the FLUKA-based (IFasso et al.ll2005l ) Monte Carlo simulation of 
the CREAM-II instrument; (open squares) from calculations usi ng the semiempirica l total 
cross-section formulae for nucleus- nucleus reactions developed by lSivher et al.l (119931 ) . 
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Fig. 5. — Survival probabilities of nuclei in a residual atmospheric overburden of 3.9 g cm 2 . 
For each elemental species, thr ee values of the s ame correction are shown, obtained respec- 
tively from the FLUKA-based jFasso et alihood ) Monte Carlo simulation of the CREAM-II 
instrument (filled circles), and from calculations using two differe nt parameteri z ation s of the 
total fra gmentation cross sections for nucleus-nucleus reactions: ISivher et al.l (119931 ) (open 
squares) ; iKox et al.l (119871 ) (open triangles) . 
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Fig. 6. — Differential intensity as a function of the kinetic energy per particle for cosmic-ray 
nuclei C, O, Ne, Mg, Si and Fe, respectively. CREA M-II data points (filled circles) are com- 
pared with previou s observations from : HEAP 3 C2 (lEngelman et al.lll990l ). trian gles; CRN 



(Mil 



ler et al.lll99lh . squares; ATIC-2 flPanov et al.ll2006f ). open circles; TRACER flAve et al. 



20081 ). stars. The solid line represents a power-law fit to the CREAM-II data. The fitted 



values of the spectral indices are listed in Table [2] 
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Fig. 7. — Energy spectra (per nucleon) for the elements C, O, Ne, Mg, Si and Fe, respec- 
tively. The differential intensities are multiplied by E 2 5 . C REAM-II results (filled circles) 
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are compar ed with previous observations by: HEAP 3 C2 (Engelman et al.lll990l ). trian- 
gles; CRN ' 



Muller et al.lll99lf) . squares; ATIC-2 flPanov et al. 



20061 ). open circles; TRACER 
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^), stars. The solid line represents a power-law fit to the CREAM-II data. 



-30- 



> 



u 



10 



X 

SIS 



10': 



_l I I 



lL 



10 



10 z 10 J 
Kinetic energy per nucleon (GeV/n) 



Fig. 8. — Measurements of the nitrogen energy spectrum (multiplied by E 2 5 ) . CREAM-II re- 



sults (filled circles) are com pared with previous observations by HE AO 3 C2 (jEngelman et al. 



1990|) (triangles) and CRN (jSwordy et al.lll990l ) (squares). The error bars in CREAM-II data 
represent only the statistical errors, while the gray bands show the systematic uncertainties 
in the differential intensity. The solid line represents a power-law fit to the combined data 
at energies > 200 GeV/n. The fitted spectral index is 2.61 ± 0.10 (x 2 /ndf = 0.4/2). 
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Fig. 9. — Measurements of the relative abundance ratios of C, N, Ne, Mg, Si and Fe, 
respectively, to O as a function of energy. CREAM data from the first (open crosses) 
and second (filled circ les) flight are compare d with previous obs ervations by HE AO 3 C2 
(jEngelman et al.lll990l ) (triangles) and CRN (jSwordy et al.lll990l ) (squares). The arrow in 
the N/O plot is an upper limit set by CRN at 1.1 TeV/n. 
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Fig. 10. — The fitted spectral indices (Table |2J) from CREAM-II data (filled circles) are 
compared with the spectral indic es of a power-law fit to: (stars) the combined CRN and 
TRACER data above 20 G eV/n ([Ave et al.l 120081): ( open squares) a compilation of direct 
and indirect measurements (IWiebel-Sooth et al.l 1 1 9981 ) . 



